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基于 遗传 算法 的 B 样 条 曲线 拟 合 改 进 算法 


高 茂 庭 ， 汉 莉 
(上 海 海事 大 学 信息 工程 学 院 ， 上 海 201306) 


摘 要 : B 样 条 曲线 拟 合 应 用 于 绘制 离散 数据 点 的 变化 趋势 ， 一 般 采 用 数据 逼近 或 者 迭代 的 方法 得 到 ， 是 图 像 处 理 和 
逆向 工程 中 的 重要 内 容 。 针 对 待 拟 合 曲线 存在 多 峰值 、 炎 点、 间断 等 问题 ， ws 
算法 。 首 先 利 用 惩罚 函数 将 带 约 束 的 曲线 优化 问题 转换 为 无 约束 问题 ， 然 后 利用 改进 的 遗传 算法 来 选择 合适 的 适应 度 
函数 ， 再 结合 模拟 退火 算法 自 适应 调整 节点 的 数量 和 位 置 ， 在 寻 优 的 过 程 中 找到 最 优 的 节点 向 量 ， 持 续 和 迭代 直到 产生 
最 终 的 优良 重建 曲线 为 止 。 实 验 结果 表明 ， 算 法 有 效 地 提高 了 精度 并 加 快 了 收敛 速度 
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Improved B-spline curve fitting algorithm based on genetic algorithm 


Gao Maoting, Feng Li 
(College of Information Engineering, Shanghai Maritime University, Shanghai 201306, China) 


Abstract: B-spline curve fitting is applied to draw the changing trend of discrete data points, which usually obtains by data 
approximation or iterative method. It plays an important part in image processing and reverse engineering. Aiming at the 
situations where multi peak, cuspidal point or discontinuity exists in the curve to fit, this paper proposed a B-spline curve fitting 
algorithm based on genetic algorithm. Firstly it used the penalty function to transform the constrained optimization problem into 
an unconstrained problem. Then it used an improved genetic algorithm to select an adaptive fitness function, and adjusted the 
number and positions of nodes adaptively by combining the simulated annealing algorithm to find the optimal node vector. The 
iterations continued until generating the final good reconstruction curve. Experimental results show that the algorithm improves 


accuracy and speeds up convergence effectively. 
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通常 采用 分 段 的 低 次 Bezier 曲线 段 进 行 拼接 , 但 是 它 缺 乏 灵 活 


性 ， 曲 线 的 次 数 严格 依赖 于 确定 该 曲线 的 数据 点 的 个 数 ， 而 且 

线 拟 合 是 用 连续 曲线 近似 地 刻画 或 比拟 一 组 离散 点 坐标 ”局 部 性 能 很 差 ， 若 要 改变 一 段 曲 线 ， 则 需要 调整 所 有 的 控制 点 

之 间 函 数 关系 的 一 种 数据 处 理 方法 ， 已 经 广泛 地 应 用 于 逆向 工 位置。 而 B 样 条 曲线 "作为 Bezier 曲线 的 一 般 化 ， 既 保留 了 

程 和 图 像 处 理 等 相关 和 领域。 曲线 拟 合算 法 主要 分 为 两 大 类 : 一 Bezier 曲线 的 优点 ， 又 增强 了 可 局 部 修改 的 特性 ， 可 容易 地 解 
类 是 参数 拟 合 ， 将 因 变 量 和 自 变量 通过 参数 方程 表示 出 来 ， 常 决 连续 性 的 拼接 问题 ， 实 现 曲 线 的 重建 。 

用 的 最 小 二 乘法 [1] 是 以 误差 的 平方 和 最 小 为 准则 , 根据 测量 数 节点 向 量 ” 和 数据 参数 ”是 决定 拟 合 效果 好 坏 的 重要 因素 。 

据 估 计 线 性 模型 中 未 知 参数 的 一 种 基本 参数 估计 方法 ， 在 理论 ”针对 了 B 样 条 曲线 拟 合 中 的 多 峰值 、 非 线性 问题 ， 以 及 待 拟 合 的 


研究 和 工程 应 用 中 都 具有 重要 的 作用 ， 同 时 又 是 许多 其 他 更 复 线 有 尖 点 、 间 断 的 情况 ， 提 出 基于 改进 的 遗传 算法 的 B 样 条 
杂 方 法 的 基础 ， 另 一 类 是 非 参 数 拟 合 ， 其 直接 找 出 因 变 量 和 自 线 拟 合算 法 ， 依 据 最 小 二 乘 模型 控制 误差 ， 并 通过 遗传 算法 
变量 之 间 的 函数 关系 ， 一 般 采 用 插值 法 。 常 用 的 有 牛顿 插值 多 ”选择 最 优 的 节点 向 量 。 在 遗传 优化 中 ， 结 合 模拟 退火 算法 "和 
项 式 法 、 差 商 与 牛顿 插值 法 、 分 段 低 次 插值 法 等 。 改进 的 遗传 算 子 克服 早熟 、 收 敛 慢 和 精度 不 高 问题 。 

Bezier 曲线 是 工程 上 较为 常用 的 一 种 拟 合 曲 线 ， 具 有 对 称 Sp 

庙 点 性 、 凸 包 性 、 几 何不 变性 、 变 差 缩 减 性 等 很 多 优良 性 1 相关 研究 
且 针 对 多 项 式 基 函 数 的 整体 性 效果 较 好 。 在 实际 应 用 中 ， B 样 条 理论 


eal 


3 

全 
J 
Me 


过 


早 由 Schoenberg 在 1946 年 提出 ,之 后 Boor、 


泾 
沽 


收 稿 日 期 : 2018-03-12; 修 回 日 期 : 2018-05-07 ”基金 项 目 : 国家 自然 科学 基金 资助 项 目 (61703267) 
作者 简介 : 高 茂 庭 〈1963-)， 男 ， 江 西 九 江 人 ， 教 授 ， 博 士 ， ee 信息 处 理 、 数 据 库 与 信息 系统 (mtgao@163.com); 冯 痢 (1994-)， 女 ， 
江苏 南通 人 ， 硕 士 研究 生 ， 主 要 研究 方向 为 数据 挖掘 .计算 机 图 形 学 


201806.00131v1 


chinaXiv: 


录用 稿 


Cox 在 1972 年 分 别 独立 定义 出 B 样 条 的 标准 算法 ,接着 Gordon 


和 Riesenfeld 把 B 样 条 运 | 
方法 ， 将 B 样 条 基 函 数 取 


于 描述 形状 ， 进 而 就 提出 了 B 样 条 
尺 Bernstain 基 函 数 ， 生 成 了 样 条 上 
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2 _B 样 条 拟 合 问题 描述 及 模型 建立 


2.1 最 小 二 乘 模型 


线 。 给 定 一 组 非 递减 序列 U: wo < Ww <… < wuntxt1i， UV 代表 节 
对 于 B 样 条 曲线 拟 合 问题 , 国内 外 的 学 者 已 经 进行 了 很 多 点 向 量 ， 那 么 B 样 条 函数 的 基 函 数 如 式 (1) (2) 所 示 。 
相应 的 研究 ， 同 时 提出 了 各 种 不 同 的 方法 。 比 如 Piegl 等 人 " 提 No = 位 0 us GD 
出 点 - 切 约束 限制 的 B 样 条 曲线 插值 算法 ; Lin 等 提出 依据 切线 t—u Urei—t 
， Nix(D) = Np (OD) + EE Neyi s(t) (2) 
的 曲率 条 件 的 三 次 B 样 条 曲线 几何 算法 ， 肖 秩 军 等 人 "提出 采 i i 
用 和 迭代 最 近 点 来 优化 修正 数据 点 的 参数 的 方法 ， 李 继 云 等 人 ” ” 式 (2) 中: kk 表示 指 向 曲线 的 阶 数 。B 样 条 的 参数 方程 为 
提出 在 给 定 的 误差 下 ,基于 最 少数 据点 的 B 样 条 曲线 拟 合算 法 ; . 
的 p(D = 》d Nu(Dss G) 
段 振 云 等 人 “结合 曲线 插值 与 B 样 条 曲线 拟 合 ， 增 加 型 值 点 以 er 
保证 曲线 局 部 性 的 算法 。 除 了 传统 的 由 方程 推导 或 者 几何 方法 式 (3) 中 : qi 表示 B 样 条 曲线 的 控制 点 坐标 向 量 ,如果 数据 点 集 
来 求解 曲线 拟 合 外 ， 也 出 现 了 多 种 结合 人 工 智能 算法 来 进行 优 。 P= fpj:j = 12,…,7] 就 在 这 条 曲线 上 ， 那 么 点 py 满足 
化 的 算法 。Yoshimoto 等 人 利用 遗传 算法 的 实数 编码 求解 无 约 了 
A i i » p(t) = Ydi Nox) (4) 
束 的 拟 合算 法 ， 孙 越 流 等 人 [9 提出 将 数据 参数 和 节点 向 量 建 3 二 
关系 ， 提 高 遗传 算法 适应 度 ， 优 化 种 群 进化 的 方法 ， 徐 善 健 等 将 式 (4) 改写 成 矩阵 形式 得 到 
人 "提出 针对 不 连续 的 曲线 ， 结 合 混 沌 蚂蚁 群 优化 算法 ， 调 整 P=N*xD (5) 
自由 节点 位 置 的 思想 ，Galvez 等 人 "在 遗传 算法 拟 合 的 基础 上 ” 式 (5) 中: P 表 示 r 个 数据 的 r x 3 的 矩阵 ，N 表 示 r xm 的 B 样 条 
运用 粒子 群 优化 动态 确定 向 量 长 度 和 节点 个 数 来 优化 曲线 等 。 基 函 数 系数 矩阵 ; d 表 示 含 有 n x 3 个 控制 点 的 控制 点 矩阵 。 当 
上 述 算法 在 一 定 程度 上 都 进行 了 优化 ， 但 精度 、 收 全 速度 和 遗 。 kk <n<r 时 , 式 (5) 可 以 用 最 小 二 乘 求 出 近似 解 , 如 式 (6) 所 示 。 
传 算法 的 “早熟 ”问题 还 需要 进一步 地 改进 。 这 些 方法 大 都 以 D=(V NIN7P (6) 
算法 的 拟 合 精度 作为 优化 目标 ， 又 综合 考虑 收敛 速度 和 人 工 智 由 此 得 到 曲线 的 误差 为 
能 算法 本 身 局 限 性 ， 从 而 将 所 求 问题 最 终 转换 成 多 约束 的 非 线 所 
E=) 00) -py 07) 
性 优化 问题 。 en 
拟 合 曲线 的 逼近 效果 很 大 程度 上 取决 于 节点 向 量 的 选择 ， E=D"(I -NN )D (8) 
合适 的 节点 向 量 对 于 曲线 的 逼近 有 着 关键 的 作用 。 本 文 利用 遗 从 式 (7) (8) 中 明显 可 以 看 出 误差 的 大 小 与 参数 的 数值 直 
传 算 法 确定 节点 向 量 ， 由 最 小 二 乘法 将 数据 参数 和 节点 向 量 联 。” 接 相关 。 
系 起 来 ， 控 制 曲线 拟 合 精度 ， 最 终 完成 曲线 的 重建 。 针 对 曲线 。 2.2 问题 提出 及 基本 解决 思路 
拟 合 中 局 部 收敛 速度 过 慢 、 容 易 出 现 早熟 现象 的 问题 ， 重 新 定 事先 确定 好 数据 参数 或 节点 向 量 的 方法 并 不 是 最 佳 的 拟 合 
义 适 应 度 函 数 ， 再 结合 模拟 退火 算 子 ， 最 后 改变 遗传 算法 的 交 方法， 选择 把 这 两 个 量 作为 变量 ， 才 是 寻求 最 好 的 拟 合 效果 的 
又 3、 变异 算 子 来 增加 种 群 的 多 样 性 ， 保 证 最 后 的 全 局 收敛 。 办 法 。 考 虑 到 两 者 之 间 的 联系 ， 根 据 最 小 二 乘 模型 来 动态 设置 
线 拟 合算 法 的 思路 流程 如 图 1 所 示 。 节点 向 量 ， 改 变 B 样 条 曲线 的 大 小 和 逼近 效果 。 由 于 所 设置 的 


曲线 拟 合 思路 流程 


节点 是 不 断 变化 的 ， 再 经 过 遗传 算法 的 迭代 之 后 并 不 能 保证 
(0) 中 的 N7TN 保 持 正定 ， 并 且 最 小 二 乘法 求 出 来 的 是 极 值 ， 不 能 
确保 一 定 是 最 值 ， 而 且 采 用 基于 函数 的 方法 求 导 也 不 能 保证 每 
个 函数 都 存在 导数 ， 所 以 选择 损失 函数 ”的 概念 来 改进 传统 的 
最 小 二 乘 解 法 。 

根据 B 样 条 曲线 的 参数 函数 ， 假 使 控制 顶点 {Pi},i = 
0,1,…,n 和 节点 向 量 {wy},j = 0,1,…,n 十 p 十 1 都 是 未 知 的 变量 ， 
那么 B 样 条 曲线 逼近 问题 就 可 以 转换 成 一 个 非 线性 的 最 小 值 优 
化 问题 号， 如 式 (9) (10) 所 示 。 


Q= mn (Pha -PP 


i=0 
m n 全 
Q=min| >》|a- > sc 
i=0 j=0 


传统 的 方法 来 求解 上 述 的 最 小 值 优化 问题 比较 困难 ， 


(9) 


(10) 


Hr 
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从 而 对 传统 方法 进行 改进 ， 通 过 遗传 算法 结合 惩罚 函数 ， 将 有  ” 是 均匀 参数 化 、 累 加 强 长 参数 法 、 向 心弦 长 参数 化 法 等 。 相 较 
约束 的 优化 问题 转换 成 无 约束 的 问题 。 在 遗传 算法 优化 中 ， 选 。 ”于 其 他 方法 ， 向 心弦 长 参数 化 "可 以 很 好 地 体现 数据 点 的 分 布 
择 一 个 合适 的 适应 度 函 数 ， 改 变 三 个 基本 算 子 ， 加 上 模拟 退火 特点， 并 且 当 数据 出 现 突然 变化 时 ， 拟 合 效果 比较 准确 。 根 据 
算 子 ， 得 到 节点 向 量 ， 接 着 带 入 最 小 二 乘 模型 中 得 到 最 优 控制 。 参数 t: 确 定数 据点 之 间 的 位 置 关 系 ， 公 式 如 式 (14) 所 示 。 


Psd 


点 ， 最 终生 成 最 优 的 曲线 拟 合 。 ty=0 
网 站 |asa — dill (14) 
3 ”改进 的 B 样 条 曲线 拟 合算 法 0 2 Na al 
针对 曲线 拟 合 对 于 光滑 、 尖 点 和 非 连 续 问题 的 处 理 不 足 ， 3.3 ”节点 向 量 的 设置 
改进 遗传 算法 的 算 子 ， 提 出 改进 的 拟 合算 法 。 算 法 的 处 理 步 又 到 目前 为 止 ， 对 于 节点 向 量 的 设置 已 经 提出 了 很 多 方法 ， 
如 下 : 主要 有 均匀 节点 向 量 . 平 均 节 点 向 量 和 皮 格 尔 逼 近 节 点 向 量 等 。 
(a) 输 入 待 拟 合 的 数据 di; 第 一 种 方法 完全 没有 考虑 数据 点 的 分 布 情况 ， 效 果 较 差 ， 第 二 
(pb) 数 据点 的 参数 化 设置 ， 以 及 相应 的 适应 度 函 数 的 建立 ; 种 虽然 考虑 了 数据 点 ， 但 是 平均 选择 过 于 简单 ， 考 虑 的 情况 较 
(c) 遗传 算法 的 参数 初始 化 , 包括 种 群 大 小 popsize、 交叉 概 少 ; 第 三 种 皮 格 尔 提 出 的 逼近 方法 ， 解 决 数 量 较 少 的 问题 还 可 
率 D-、 变 异 概率 pn 、 最 大 迭代 次 数 ; 以 ， 但 一 旦 数据 量 增 大 ， 计 算 量 增加 ， 算 法 效率 不 高 。 针 对 这 
(gd) 根据 格雷 码 的 编码 方式 初始 化 种 群 8(b); 种 情况 ， 本 文 提 出 新 的 选择 方式 ， 有 具体 步骤 如 下 : 
(e) 由 最 小 二 乘法 得 出 曲线 控制 点 ,并 计算 每 一 个 个 体 的 适 a) 对 端点 的 控制 点 进行 插值 。 一 个 非 递减 的 节点 向 量 vo < 
应 度 函 数值 ; U1 三 … < Untp+t1, 端点 节点 的 重复 度 设 置 成 B 样 条 的 阶 数 为 
(f) 迭代 计算 , 最 终 输 出 最 优 的 节点 向 量 和 控制 点 , 重建 p+1:uw=W = =Up=0, Un = Un2 = = Unp+ti = 1; 
线 ， 否 则 转 到 (8g) ; b) 把 首尾 位 置 的 数据 点 作为 重建 的 B 样 条 曲线 的 首尾 控制 
(g) 对 6(t) 由 上 文 的 选择 算 子 执行 选择 操作 得 61(t); 点 ， 即 do = Pu dm = PP; 
(hb) 对 61(t) 由 上 文 的 交叉 算 子 执行 交叉 操作 得 62(t); c) 除 了 首尾 外 剩 下 了 n 十 p 个 内 部 的 节点 , 这 些 对 于 拟 合 之 
(i) 对 52(t) 由 上 文 的 变异 算 子 执行 变异 操作 得 63(t)， 然 后 前 都 是 未 知 数 ， 运 用 遗传 算法 和 迭代 进行 自 适 应 设置 ， 最 终生 成 
转 (e) 。 最 优 节点 向 量 为 止 。 
依据 上 述 的 算法 步骤 ， 重 点 改进 最 小 二 乘法 和 遗传 算 子 得 ”3.4 遗传 算法 原理 及 相应 改进 
到 最 优 的 节点 向 量 ， 接 下 来 详细 介绍 具体 过 程 。 3.4.1 遗传 算法 原理 
3.1 传统 最 小 二 乘 解法 的 改进 遗传 算法 ”是 借鉴 大 自然 的 “ 适 者 生存 ”优胜 劣 汰 ”的 思 
在 最 小 二 乘法 中 ， 可 以 将 最 终 的 误差 表达 式 理解 成 损失 函 想 ， 最 美国 人 Holland 教授 在 他 的 《自然 界 和 人 工 系统 的 
数 ， 需 要 使 用 梯度 向 量 和 雅 可 比 向 量 (Jacobian matrix) 表示。 适应 性 》 中 提出 ， 是 一 种 全 局 优化 的 自 适应 概率 搜索 算法 。 遗 
定义 雅 可 比 矩 阵 是 以 误差 对 参数 的 偏 导数 为 元 素 ， 如 式 (11) 所 传 算法 以 个 体 的 适应 度 函 数值 为 依据 ， 再 在 进化 过 程 中 进行 选 
示 。 择 、 交 又 和 变异 操作 ， 完 成 寻找 问题 的 最 优 解 的 过 程 。 它 提供 
de 了 一 种 求解 复杂 系统 优化 问题 的 通用 框架 ， 在 生物 技术 、 化 学 
If W) = FE = 2, mB) = 12,,) GD ae ge 
7 工程 、 计 算 机 辅助 设计 、 医 学 工程 等 方面 都 有 着 重要 的 应 用 。 


式 (11) 中 : m 表 示 数 据点 的 数量 ，n 表 示 神 经 网 络 的 参数 数量 ， ”3.4.2 算 子 改进 以 及 参数 设置 

那么 雅 可 比 和 矩阵 就 是 m x n 阶 矩阵 。 损 失 函 数 的 梯度 向 量 就 如 简单 遗传 算法 主要 包括 初始 种 群 的 选取 、 染 色 体 的 编码 、 

式 (12) 所 示 。 适应 度 函 数 的 选取 、 遗 传 操 作 的 设计 和 算法 控制 参数 的 设 定 五 

Vf =2]"*e (12) 个 部 分 。 接 下 来 进行 具体 描述 。 
式 (12) 中 : e 表 示 所 有 的 误差 的 向 量 值 。 对 于 权 值 函数 的 更 新 和 a) 初始 种 群 的 选取 。 
优化 如 式 (13) 所 示 。 常用 的 方式 都 是 随机 生成 初始 种 群 ， 但 考虑 到 初始 种 群 的 
Wir1 = wi — UF *Jit il) ! * (2J{ * ei) (13) 设置 对 于 结果 和 算法 效率 都 是 相当 重要 的 ， 而 随机 会 使 得 个 体 

式 (13) 中 : s 表 示 衰 减 因 子 ， 用 来 保证 矩阵 的 正定 。 分 布 在 可 行 解 的 空间 是 不 均匀 的 ， 因 此 选择 均匀 设计 。 首 先 将 

3.2 ”数据 点 参数 化 的 设置 解 空 间 分 成 Z 个 区 间 ; 然后 用 均匀 数组 来 选择 P 个 染色 体 ; 最 后 
以 最 小 二 乘 为 拟 合 的 基础 ， 那 么 影响 曲线 形状 的 因素 有 两 从 Z xP 中 选择 适应 度 值 较 高 的 Q 个 个 体 作 为 初始 种 群 ， 从 而 确 

个 : a) 数 据点 参数 化 的 分 布 ，b) 节点 向 量 的 选取 。 不 同 的 数据 。”” 保 个 体 在 解 空间 中 尽量 分 散 。 

点 的 参数 和 节点 向 量 对 曲线 形状 有 着 关键 的 作用 ， 不 合理 的 选 bp) 染色 体 的 编码 。 

取 会 导致 丙 近 效果 很 差 。 遗传 算法 常用 的 的 染色 体 编码 方式 是 二 进 制 或 十 进 制 编码 。 


到 目前 为 止 ， 数据 点 参数 化 的 确定 方法 很 多 ， 比 较 常用 的 。” ”在 进行 数值 优化 时 ， 上 有 具有 稳定 性 高 、 种 群 多样 性 大 的 优点 ;但 
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其 随机 性 使 得 局 部 搜索 能 力 较 差 。 比 如 对 于 一 些 高 精度 的 问题 ， 
当 解 接近 于 最 优 解 后 , 由 于 其 变异 后 表现 型 变化 很 大 ,不 连续 ， 
所 以 会 远离 最 优 解 ， 达 不 到 稳定 。 由 此 本 文选 取 格雷 码 ， 其 连 
续 两 个 整数 所 对 应 的 编码 值 之 间 仅 仅 只 有 一 个 码 位 是 不 同 的 。 
进 制 编码 转 格雷 码 的 转换 公式 为 
gm = bm (15) 
gi=bi+1@bi(i= 
c) 适 应 度 函 数 的 选取 。 
适应 度 函 数 是 遗传 算法 中 其 他 关键 步骤 的 依赖 基础 ， 选 择 
一 个 合适 的 适应 度 函数 是 遗传 算法 优化 的 关键 。 依 据 前 文 提 出 
的 最 小 二 乘 模型 ， 要 得 到 优良 的 曲线 拟 合 结果 ， 就 要 控制 误差 
值 ， 并 且 要 减少 控制 点 个 数 。 为 此 ， 可 以 设计 适应 度 函 数 如 式 
(17) 所 示 。 


Tt 


m—1,m—2,.…,1) (16) 


fitness 5 CE (17) 
式 (17) 中 : 9 表示 控制 因子 ,表示 修正 因子 。 再 考虑 遗传 算法 


的 迭代 并 不 能 保证 节点 的 顺序 ， 而 无 序 的 节点 是 无 法 生成 B 样 
条 曲线 的 ， 可 以 设计 一 个 函数 体 。 具 体 如 下 : (a) 先 对 所 有 的 个 
体 按照 大 小 顺序 排序 ; (b) 然后 对 以 上 的 数据 点 采取 最 小 二 乘法 
求解 控制 点 ; (c) 然后 根据 适应 度 函数 确定 个 体 的 适应 度 值 。 通 
过 函数 体 的 思想 ， 将 每 次 迭代 的 结果 放 在 一 起 ， 保 证 节点 向 量 
的 有 序 性 。 

d) 遗传 算 子 的 设计 。 

选择 算 子 : 主要 体现 适 者 生存 的 原理 。 由 适应 度 函 数值 大 
小 确定 ， 保 留 适应 度 值 较 大 的 优秀 个 体 ， 丢 弃 适 应 度 值 较 小 的 
个 体 。 简单 遗传 算法 中 , 轮 盘 赌 选 择 应 用 较 广 , 但 随机 性 太 大 ， 
为 此 ， 采 用 基于 竞争 指数 的 模拟 退火 选择 算 子 ”。 采 用 竞争 指 
数 ] 综 合 考虑 个 体 的 适应 度 值 fitness， 简 单 记 为 F 和 编码 的 差异 
度 值 C， 如 式 (18) 所 示 。 
J=mrFr+m2C (18) 


式 (18) 中 : ms 和 ms 表示 比例 参数 ， 和 为 1。 将 个 体 按照 竞争 


XY 例 作 期 刊 ， 


条 曲线 和 


天 = 一 > @ logm (20) 


其 中 : pi = 锅 ， 如 果糖 数值 较 大 表示 种 群 的 种 类 较 多 ， 进 化 能 
力 较 强 ， 交 叉 概 率 就 需要 取 较 小 数值 ， 反 之 取 较 大 数值 。 最 后 
交叉 概率 如 式 (21) 所 示 。 

由 于 初始 的 数值 都 是 按照 递增 的 顺序 , 再 考虑 到 编码 方式 ， 
所 以 采用 算术 交叉 。 设 两 个 个 体 分 别 为 如 和 xz， 那 么 交叉 公式 
为 


Xi = O1X1 + (WaX2 (22) 
X2 = WiXz + WaX1 (23) 
式 (22) (23) 中 : w 表 示 (0,1) 之 间 的 随机 数 ; wi 和 ws 和 为 1。 


变异 算 子 : 变异 概率 是 直接 影响 遗传 算法 优化 结果 的 重要 
因子 。 概 率 取 值 不 能 大 范围 地 破坏 优良 基因 ， 并 且 最 好 可 以 生 
成 优良 个 体 。 在 进行 变异 操作 时 ， 适 应 度 值 低 于 平均 值 的 个 体 
被 定义 为 劣质 个 体 ， 需 要 以 较 高 的 概率 进行 变异 。 根 据 上 文 的 
思路 ， 可 以 加 入 早熟 判断 标准 ， 考 虑 当前 个 体 的 最 大 适应 度 值 
与 超过 平均 适应 度 个 体 的 平均 适应 度 值 之 差 9。 变 异 概 率 定义 
为 则 如 式 (24) 所 示 。 


1 


1 + ee0) C4) 


pm 一 


式 (24) 中 : k 是 控制 因子 ， 大 于 0， 变 异 概率 与 6 同步 变化 。 满 
足 变异 条 件 时 ， 在 每 个 基因 疡 中 产生 一 个 随机 数 s， 如 果 这 个 


随机 数 小 于 变异 概率 ， 执 行 单 点 变异 ， 如 式 (25) 所 示 。 
U*jz 
ji= Jitu* (i ji) (25) 
Oit+w 
式 (25) 中 : wu 表示 (0, 1) 中 的 随机 数 ， 秉 持 了 最 优选 择 和 最 差 


抛弃 原则 。 

如 果 持 续 按 照 上 述 步骤 执行 ， 尽 管 结果 可 以 尽快 向 最 优 解 
靠近 ， 但 是 个 体 之 间 的 相似 度 会 越 来 越 高 ， 种 群 的 多 样 性 也 就 
逐步 降低 ， 很 容易 出 现 局 部 收敛 的 现象 。 为 了 解决 这 个 问题 ， 


指数 的 大 小 排序 ， 选 择 的 先后 顺序 与 位 置 先后 一 致 ， 合 理 保 持 
群体 优良 性 。 
交叉 算 子 : 好 的 交叉 操作 可 以 使 个 体 之 间 的 优良 基因 进行 
信息 交换 。 在 进行 交叉 操作 时 ， 适 应 度 函数 值 高 于 平均 适应 度 
值 的 个 体 被 定义 为 优良 个 体 ， 需 要 以 较 高 的 概率 交叉 ， 以 保证 
保留 优良 基因 。 但 如 果 仅 依靠 平均 适应 度 值 ， 易 陷入 局 部 优 
化 ， 因 此 考虑 加 入 方差 和 焙 来 设计 概率 。 
妥 设 种 群 的 规模 大 小 是 W， 个 体 适 应 度 值 是 卢 ， 平 均值 是 
_f， 方 差 表 达 如 式 (19) 所 示 。 


G1= 2 [Gf /N] (19) 


方差 代表 群体 中 个 体 的 分 布 情况 ， 当 方差 较 大 即 个 体 较 为 
分 散 时 ， 减 小 交叉 概率 ， 反 之 增加 。 

炉 是 将 问题 的 可 行 解 分 为 4; 个 范围 ， 那 么 种 群落 在 4; 的 个 
体 数目 为 |4i|， 和 如 式 (20) 所 示 。 


结合 上 文 的 早熟 判断 指标 ， 达 到 指标 时 ， 过 滤 到 相似 个 体 。 具 

体 做 法 如 下 : 
a) 将 种 群 个 体 按照 适应 度 值 进行 排序 ; 

) 根 据 编码 方式 计算 相似 度 ， 确 定 一 个 B 值 ; 

c) 把 相似 值 结果 与 B 进 行 比较 ， 删 除 相 似 的 个 体 。 

完成 过 滤 操作 后 ， 需 要 有 新 的 个 体 来 补充 ， 不 能 用 直接 随 

E 成 的 方式 ， 不 利于 全 局 搜索 ， 可 以 将 原本 中 适应 度 比 较 高 

的 个 体 进行 变异 ， 增 加 优良 个 体 ， 与 其 他 原本 的 个 体 也 更 易 产 

生 有 优良 的 子 代 。 


4 ”实验 结果 与 分 析 


为 了 检验 本 文 提出 的 基于 遗传 算法 的 B 样 条 曲线 拟 合算 法 
的 实际 应 用 效果 ， 针 对 具有 平滑 、 不 连续 和 有 人 尖 点 特点 的 测试 
函数 ， 选 择 出 合适 的 节点 ， 然 后 得 出 拟 合 结果 ， 并 与 传统 简单 


遗传 算法 、 结 合 蚁 群 算法 优化 的 拟 合 方法 进行 比较 。 
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4.1 实验 环境 

实验 硬件 环境 是 CPU Intel® CoreCTM 2 Duo, 内 存 为 4GB， 
操作 系统 环境 是 Windows7， 编 程 平 台 和 实验 环境 为 MATLAB 
8.0。 
4.2 ”测试 函数 的 选择 
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一 一 迭代 次 数 和 节点 数 关系 表示 


考虑 到 具有 平滑 、 不 连续 和 有 人 尖 点 特点 的 测试 函数 ， 可 以 
测试 出 曲线 拟 合算 法 的 多 适应 性 。 采 用 文献 [11] 中 的 三 个 函数 
作为 测试 函数 ， 这 也 是 较 多 B 样 条 曲线 拟 合 选择 的 函数 ， 它 们 
的 表达 式 和 定义 域 如 表 1 所 示 。 1 人 is0 200 
表 1 测试 函数 送 代 次 数 
函数 表达 式 定义 域 nk 4 ”测试 函数 户 Co 的 迭代 次 数 与 节点 数 的 关系 
ee 0 通过 对 图 (2) 一 (4) 分 析 可 得 ， 卢 (的 最 优 节点 数目 为 4 
， 三 (x) 的 最 优 节点 数目 为 8， 包 (x) 的 最 优 节点 数目 为 5。 
[Ae | 001 + G03) Te 4.3.2 测试 函数 拟 合 效果 
DIS TC 0.65) 2 由 求 出 的 最 优 节点 ， 进 行 曲线 拟 合 ， 三 个 测试 函数 的 拟 合 
R(x) i xe[od] 1~17 效果 如 图 (5) 一 (7) 所 示 。 


4.3 实验 结果 及 分 析 100 
4.3.1 最 优 节点 数目 的 选取 


80 

通过 对 测试 函数 请 (x)、 户 (x) 和 太 (x) 的 节点 数目 随 着 迭代 | 
次 数 的 增加 而 发 生 的 演变 趋势 的 分 析 ， 从 而 可 以 得 到 最 优 的 节 
点 数目 ， 重 建 出 函数 曲线 的 拟 合 图 。 三 个 测试 函数 的 迭代 次 数 0 
和 节点 数目 的 函数 关系 如 图 CD) 一 (4) 所 示 。 ol 

20 T T T 0 
迭代 次 数 和 节点 数 关系 表示 0 2 0 6 oe 1 
图 5 测试 函数 有 (x) 的 拟 合 效 果 


拟 合 曲线 
0 上 1 LL | 
0 50 100 150 200 
迭 代 次 数 
图 2 测试 函数 (x) 的 迭代 次 数 与 节点 数 的 关系 
0 i . i i 
25 0 02 0.4 0.6 0.8 1 
一 一 一 迭代 次 数 和 节点 数 关系 表示 


图 6 测试 函数 户 (的 拟 合 效果 


5 
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迭代 次 数 本 
图 3 ”测试 函数 户 (Ca 的 从 代 次 数 与 节点 数 的 关系 -20 


0 0.2 0.4 0.6 0.8 


图 7 测试 函数 户 (x) 的 拟 合 效果 
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(5) 一 (7) 分 别 表示 运用 本 文 提出 的 方法 对 测试 函数 的 拟 
合 效果 , 图 (5) 是 带 有 平滑 部 分 的 曲线 拟 合 效果 , 从 中 可 以 看 出 ， 
在 平滑 位 置 的 拟 合 效果 较 好 ， 虽 然 部 分 区 域 存 在 一 定 的 偏差 ， 
但 都 在 可 接受 范围 内 图 (6) 是 不 连续 的 曲线 拟 合 情 况 ， 由 图 可 
以 看 出 在 间断 点 处 的 处 理 很 好 ， 接 近 原 来 的 曲线 趋势 ; 图 (7) 是 
尖 点 的 曲线 ， 在 拟 合 的 图 中 可 以 明显 有 极 值 的 情况 ， 而 且 基 
本 展示 出 其 原本 的 特征 。 通 过 以 上 三 个 测试 函数 的 拟 合 效果 ， 
表明 基于 改进 遗传 算法 的 B 样 条 曲线 拟 合算 法 对 于 平滑 、 不 连 
续 和 带 尖 点 三 种 情况 处 理 良 好 ， 拟 合 效 果 与 实际 曲线 趋势 几乎 
一 致 ， 由 此 说 明 算 法 性 能 优良 。 
4.3.3 与 其 他 拟 合算 法 的 比较 

为 了 进一步 评估 本 文 提出 的 算法 的 合理 性 ， 对 比分 析 三 种 
算法 按 式 (8) 计算 结果 。 表 2 给 出 了 本 文 算法 与 文献 [1] 、 文 献 
[11] 的 拟 合 误差 。 


| 
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表 2 测试 函数 的 误差 比较 
测试 函数 fi(x) 的 数据 点 。 ”fz() 的 数据 点 “fs(9) 的 数据 点 
误差 误差 误差 
本 文 方法 1.075e-06 4.062e-06 1.773e-02 
文献 []] 的 方法 5.967e+00 9.327e-01 7.762e+01 
文献 [11] 的 方法 2.563e-04 3.743e-04 5.532e-02 
表 (2) 的 数据 可 知 ， 三 个 测试 函数 在 应 用 本 文 算法 拟 合 
时 ， 误 差 都 相对 较 小 ， 逼 近 效 果 更 好 。 
5 ”结束语 
本 文 从 遗传 算法 的 重要 算 子 进行 改进 , 结合 最 小 二 乘 模型 ， 
实验 结果 表明 ， 采 用 上 述 的 改进 遗传 算法 对 We 
更 好 ， 精 度 更 高 。 但 是 对 于 遗传 算法 要 获得 更 高 的 数值 精度 ， 


就 需要 提高 相应 的 收敛 速度 ， 而 且 数 值 的 参数 化 还 需要 加 以 提 
升 ， 最 后 遗传 算法 的 随机 性 很 大 ， 增 加 算法 的 稳定 性 将 是 以 后 
值得 继续 研究 的 问题 。 之 后 还 可 以 继续 将 目前 研究 的 二 维 曲 线 
问题 拓展 到 三 维 曲 面 上 ， 从 而 紧 跟 时 代 发 展 的 步伐 。 
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